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\  ABSTRACT 

This  paper  introduces  a  contour-based  approach  to  motion 
estimation.  It  is  based  on  first  computing  motion  at  image 
corners,  and  then  propagating  the  corner  motion  estimates  along 
the  principal  contours  in  the  image  based  on  a  local  2hD  motion 
assumption.  The  results  of  several  experiments  are  presented. 
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1.  Introduction 

The  earliest  problem  that  arises  in  the  analysis  of  time- 
varying  images  is  the  detection  of  moving  image  elements 
(edge,  regions)  and  the  computation  of  the  image  velocity 
(optic  flow)  of  those  elements.  A  variety  of  computational 
schemes  have  been  proposed  to  solve  this  problem.  In  a 
recent  survey,  Ullman  [1]  broadly  classifies  these  as  inten¬ 
sity-based  and  token-matching  schemes. 

An  important  class  of  intensity-based  schemes  takes  ad¬ 
vantage  of  the  relationship  between  the  temporal  and  spatial 
gradient  of  any  continuous  and  differentiable  image  proper¬ 
ty  which  is  invariant  to  small  changes  in  perspective.  For 
example,  if  we  assume  that  the  intensity,  I,  satisfies  these 
properties,  the  relationship 

-It  *  “Ix  *  ''ly 

can  be  used  to  determine  velocity.  Here,  is  the  temporal 

intensity  gradient,  I^  and  I  the  x  and  y  components  of  the 

X  y 

spatial  intensity  gradient,  and  u  and  v  the  x  and  y  compo¬ 
nents  of  image  velocity.  Measuring  I. ,  I  ,  I  from  an  image 

u  X  y 

sequence  establishes  a  linear  constraint  on  the  x  and  y  velo¬ 
city  components.  A  single  velocity  estimate  can  be  computed 
by  spatially  combining  the  constraints  using  e.g..  Hough 
transforms  (2] ,  least-squares  methods  [3]  or  minimization 
techniques  [4] .  All  of  these  techniques  suffer  from  certain 
disadvantages.  The  Hough-transform  and  minimization  techniques 


assume  that  image  velocity  is  uniform  over  large  parts  of 


the  image,  and  the  least-squares  method  further  assumes  that 
the  constraint  equations  determined  for  nearby  points  are 
independent  -  an  assumption  that  is  violated  by  the  spatial 
integration  required  to  compute  spatial  derivatives. 

In  this  paper  we  develop  a  contour-based  approach  to 
motion  estimation  at  a  small  set  of  image  points  at  which 
it  is  possible,  in  principle,  to  unambiguously  determine 
image  velocity.  Specifically,  corners  have  the  property 
that  their  motion  can  be  directly  computed  based  only  on 
measurements  made  at  the  corner  (in  practice,  of  course, 
one  must  examine  a  small  neighborhood  of  the  corner). 

Another  important  property  of  corners  is  that  they  can 
be  safely  regarded  as  projections  of  scene  features  whose 
general  appearance  is  invariant  to  rigid  motion  -  e.g.,  an 
image  corner  may  be  the  projection  of  the  vertex  of  a  poly¬ 
hedron,  or  of  a  curvature  discontinuity  on  the  boundary  of 
a  surface  marking.  The  second  step  involves  propagating 
these  velocity  estimates  to  a  larger  number  of  picture  points. 
This  is  based  on  the  assumption  that  the  image  motion  is 
locally  a  rigid  two-dimensional  motion.  Given  this  assump¬ 
tion,  the  velocity  at  one  point  on  a  contour  and  the  normal  com¬ 
ponent  at  a  neighboring  point  can  be  combined  to  compute  the 
actual  velocity  at  the  neighboring  point. 


2.  Estimating  motion  at  corners 


The  motion  of  a  corner  can  be  computed  in  a  variety  of 
ways.  Section  2.1  describes  an  approach  based  on  temporal 
intensity  changes  along  lines  parallel  to  the  sides  of 
the  corner.  Section  2.2  discusses  a  second  technique 
which  combines  normal  vectors  in  a  small  neighborhood  of 
the  corner  along  the  contour.  It  is  similar  to  the  velo¬ 
city  estimation  algorithm  in  Horn  and  Schunck  [4].  We 
should  also  point  out  that  Nagel  [5]  has  recently  proposed 
a  corner  velocity  estimation  algorithm  based  on  a  differ¬ 
ential  approach  (i.e.,  a  Taylor  series  expansion  of  the 
image  function  in  the  neighborhood  of  the  corner  truncated 
after  the  second-order  terms.). 

2.1  A  structural  approach 

This  subsection  presents  a  structural  approach  to  cor¬ 
ner  motion  estimation.  We  first  describe  velocity  compu¬ 
tation  for  the  case  of  translation  motion,  and  then  con¬ 
sider  translation  combined  with  rotation. 

Suppose  that  a  corner  simply  translates  from  point  Cq 
to  between  rwo  frames  t^  and  t-^  (see  Figure  1).  Let  O' 
be  a  point  on  the  bisector  of  PCqR  and  let  O' A  and  O'B  be 
lines  parallel  to  CqP  and  CqR,  respectively  at  some  unit 
distance  from  C^P  and  CqR.  Suppose  that  |0'A|  =  lO'B]  = 

1+m,  for  some  constant  m.  Finally,  assume  that  the  intensity 
inside  the  corner  is  1  and  outside  the  corner  is  0. 


Now,  at  time  t^,  the  average  intensity  along  line  seg¬ 
ments  O' A  and  O'B  is 

”  l/(l+m) 

If  Ax'  and  Ay'  are  the  components  of  the  translation  in 
the  directions  of  the  lines  O'A  and  O'B,  then 

^O'A^^l^  ~  d+Ax' )/(l+m) 

^O'B^^l^  =  d+Ay' )/(l+m) 

assuming  that  m  is  chosen  large  enough  so  that  max (Ax', Ay')  < 
m.  Finally,  Ax'  and  Ay'  can  be  computed  from 

^O'A  “  ^O'A^^l^  ~  ^O'A^^O^  ”  Ax'/d+m) 

^O'B  ~  ^O'B^^l^  ~  ^O'B^^'O^  '  Ay'/d+m) 

Once  Ax'  and  Ay'  are  computed,  the  components  of  the 

velocity  in  the  original  image  coordinate  system  can  be 
recovered  easily: 

[Axl  Fcosa  cosgl  r^'l 
AyJ  [sina  sin&J  lAy'J 

The  practical  success  of  this  technique  depends  on  our 
ability  to  compute  several  corner  parameters  accurately. 

These  parameters  are 

1.  corner  location  at  tp, 

2.  corner  shape  Cangles  a  cuid  3)  ,  and 

3.  corner  contrast  Cassumed  here  to  be  1) 

The  computation  of  these  parameters  is  discussed  in  Section  4.1. 
Next,  we  extend  the  previous  simple  analysis  to  include  rota¬ 
tion  as  well  as  translation.  We  will  treat  this  case  as  a 


translation  from  Cq  to  followed  by  a  rotation  about 
through  a  clockwise  angle  y  (see  Figure  2) .  Since  translation 
and  rotation  are  specified  by  a  total  of  three  parameters, 
we  could  extend  the  above  analysis  using  only  a  third  line 
segment  parallel  to  either  O' A  or  O'B.  Instead,  we  consider 
two  pairs  of  parallel  line  segments,  and  compute  the  displace¬ 
ments  in  the  directions  O'A  and  O'B  rather  than  directly 
computing  the  angle  y. 

Let  Ax'^,  Ay'^  be  the  translational  components  of  the 
motion  in  the  O'A  and  O'B  directions,  and  Ax'^.  and  Ay'^ 
the  corresponding  rotational  components.  Then 

Ax'^  +  Ax'j.  =  (l+m)  AIq,^  (2.1) 

Ay'^  -  Ay'j.  =  (l+m)  AI^.^  (2.2) 

From  Figure  3,  we  see  that 


^  x-  -  (2.3) 

1+Ay't  ^ 

where  5  is  the  dist^mce  between  the  parallel  line  segments 
O'A  and  CD.  Similarly 

Ay'r-  =  Ay'^  -  Ay"^  (2.4) 

1+Ax'^  6 

Also 

Ax'j.  -  Ax”^  =  (l+ra)IlQ,^(t^)-I^jj(t^)]  (2.5) 

Ay'j.  -  Ay"^  =  (l+m)  IlQ,g(tj^}-Ij.j.(t^)]  (2.6) 

Substituting  (2.5)  aixi  (2.6)  into  (2.3)  and  (2.4)  and  simplify¬ 
ing  ,  we  obtain 

Ax'  -  C,Ay'  =  C,  (2.7) 

r  I  t  1 


C2Ax't  ”  Ay'^ 


(2.8) 


B 


where 


_  tl+ml 
^1  6 

^  _  (l+ml 

^2  6 

solving  for  Ax'^  and  Ay'^  we  obtain 

Ax'  =  (Itm)  (AIq.a  -  CiAIn.p)  -  Ci  (H-C,)  ^2.95 


I+C1C2 


Ay 


,  ^  Cl-»-mHC2AJQ,^  -  ~  (2.10) 


I+C1C2 


Substituting  (2.7)  and  (2.8)  into  (2.9)  and  (2.10),  we  can  also 
compute  Ax'^  and  Ay'^,  which  gives  us  a  complete  description 
of  the  motion  of  the  corner. 

« 

2 . 2  A  least-squcures  approach 

In  this  section  we  show  hov;  simple  least-squares  algo¬ 
rithms  can  be  used  to  compute  corner  motion.  We  can  rewrite 
(1.1)  to  obtain 

=  -I^/lVlj  (2.11) 

where  |Vl|  is  the  magnitude  of  the  spatial  intensity  gradient 

at  that  point  (|Vl|  =  and  V  is  the  projection  of  the 

X  y 

velocity  vector  V  onto  the  intensity  gradient  at  that  point. 

We  will  first  consider  the  case  when  the  velocity  is  only  a 
translation,  and  then  consider  translation  with  rotation. 

If  we  assume  that  the  velocities  are  constant  in  a  small 
neighborhood  of  the  corner  along  the  contour,  then  we  can 
relate  the  problem  of  determining  the  velocity  V  at  the  cor¬ 
ner  to  that  of  determining  a  V  to  minimize 


(V-ff.  -  V  .  )^}  (2.12) 

~  1  ni 

where  ^  unit  normal  vector  of  the  con¬ 
tour  point  and  is  the  projection  of  V  onto  the  intensity 

gradient  direction  at  the  i^^  contour  point,  called  the  normal 
projection  for  short.  By  rainimizing  the  error  E  , 
we  obtain 


where 


au  +  cv  =  d 
cu  +  bv  =  e 


(2.13) 


a  =  Z  n^, 
i=l 

^  2 

b  =  Z  nt_ 
i=l  ^2 


c  =  Z  n. , -n. ^ 

i=l  12 


d  =  Z  n.  ‘V  . 

i=l  11  ni 


e  =  Z  n.^*V  . 

i=l  12  ni 


From  equation  (2.12),  we  have  the  velocity  estimation 


of  a  turning  point 


u 


bd-ce 

K  2  ' 

ab-c 

ae-cd 

T  2"  • 

ab-c 


(2.14) 


V 


Notice  that  the  solutions  for  u  and  v  are  only  meaningful  v;hen 
ab-c^,  which  is  related  to  the  variance  of  normal  directions,  is 
high.  For  a  straight  line  segment,  e.g.,  there  is  no  solution 
because  the  denominators  of  eqs.  (2-13)  are  zero. 


ab-c  = 


k  _  k  2 

i=l  i=l  i=l  ^ 


i=l 


=  kn^*kn2  -  (kn^^-n^)^ 


=  0 


i.e.,  from  a  small  element  of  a  straight  line,  the  only  infor¬ 
mation  that  one  can  obtaiji  is  the  motion  component  normal  to 
that  line.,  and  motion  along  this  line  element  cannot  be  detected. 
Corners,  however,  are  just  those  points  where  the  variance  of 
normal  directions  is  locally  maximal. 

Now,  consider  the  case  in  which  the  motion  of  the  contour 
can  be  decomposed  into  a  translation  with  velocity  Vq  ^^Ox'^^Oy^ 
at  the  corner  (xQ,yQ)  and  a  rotation  around  (Xp,yQ)  with  the  an¬ 
gular  velocity  w,  as  shown  in  Figure  4.  Note  that  since  we 
are  only  interested  in  the  motion  of  the  corner,  we  do  not  ex¬ 
plicitly  solve  for  w,  although  it  would  be  easy  to  do  so. 

Vie  obtain 

w  X  di  =  Vi  -  Vq  (2.15) 

By  rewriting  this  equation,  we  have 


V. 
— 1 


(2.16) 


where  components  of  the  displacement  vector 

from  point  (J^q/Yq)  to  point  (x^,y^)  on  the  contour.  The 

normal  component  of  is  related  to  d^,  w  and  by 

V  .  =  V. -n, 
ni  — i  1 


=  (V-  -wd.  )n.  +  (V-  +wd.  )n. ,  (2 

'Ox  ly  ix  Oy  ix  ly 

By  considering  three  points  on  the  contour,  and  can 


(2.17) 


be  simply  obtained  by  solving  linear  equations.  In  general, 
more  than  three  image  points  are  taken  and  ^Qx'^Oy  computed 
by  minimizing  the  following  square  error; 


=  T  (V,,  n.  +V-  -n.  +w(d.  n.  -d .  n.  )-V  .)^ 
^  Ox  IX  Oy  ly  ix  ly  ly  ix  ni 


We  obtain  the  following  least  squares  solution  of  • 


1 

^11 

a 

2 

^02 

a 

3 

‘^01 

a 

6 

'20 

a 

'll 

^2 

a, 

‘10 

S 

a 

(2.18) 


where 


S  ^  =  r.  n^  n?  d'’  d^ 

rpqk  ^  IX  ly  ix  ly 

=  X  n^  n? 

rp  ^  IX  ly 


and 


®1110"^2001 


and 


= 

^  ^0210”®1101 
^OllO'^lOCl 


c,  = 


H  V  . -n.  , 

^  XI  IX 


Co  = 


Z  V  . -n.  , 
i  ni  ly 


C-,  = 


E  V  .  , 
i 


^20  ^11  ^1 

®11  ^02  ^2 

^10  ®01  ^3 


3 .  Propagation  of  velocity  vectors  along  image  contours 


3. 1  The  local  constraint  and  the  propagation  formula 

Suppose  the  velocity  vectors  at  the  ends  of  a  contour 

A_A.  are  known  (see  Figure  4).  Consider  a  small  line  segment  ds 
U  K 

along  the  contour  AqA^^.  Assuming  that  the  motion  is  a  rigid  motion 
Vpg  of  Vjj,  the  motion  of  A^,  parallel  to  AqAj^  must  equal  the 
parallel  component  V^g  of  the  velocity  at  Aj^ : 


V 


OS 


IS 


(3.1a) 


or 


Vq  •  as  =  •  dS.  (3.1b) 

where  Vq  and  are  the  velocity  vectors  at  the  two  ends  of 
the  line  dS,  and  dS  is  the  unit  vector  along  dS,  the  vector 
joining  Aq  to  Aj^.  Rewriting  this  local  constraint  (eq.3.1b) 
into  component  form,  we  obtain 

Vq  •  dS  =  (V^^n  +  V^^t)  •  dS 

=  .  d^  +  V^^t  •  dS  (3.2) 

where  and  are  the  normal  component  and  the  tangential 
component  of  the  velocity  vector  respectively,  and  n  and 
t  are  the  unit  vectors  in  the  normal  and  tangent  directions 
of  the  contour  at  A^^.  From  Figure  5,  we  see  that 


n'dS  =  cosa  and  t.dS  =  cos(Tt/2-a)  -  sina  so  that 


'os 


=  V^^sinnt  +  V^^cosa 


(3.3) 


In 


Thus,  we  have  that  the  tangential  component  is 


=  CVgg  -  •  cosa)/sina  (3.4) 

where  a  is  the  amgle  between  the  unit  vector  dS  and  the 

normal  vector  n  at  the  point  Aj^.  We  also  have  Y=&-a,  where 

3  is  the  angle  between  the  x~axis  and  the  normal  vector  ii, 

and  Y  is  the  angle  between  the  x-axis  and  the  line  segment  dS. 

We  can  propagate  the  velocity  along  a  contour  using  eq.  (3.4), 

because  the  first  projection  V^g  is  known  after  the  previous 

propagation  and  the  normal  component  can  be  computed  by, 

e.g.,  the  methods  discussed  in  [3]  or  [4],  Once  :.s  con- 

*1  6 

puted,  Vj^=Vj^e-’  can  be  obtained  because 


''i  = 


.3.5) 


0=3-  arctan  V.^/V. 

It'  In 


From  eq.  (3.4)  the  new  estimate  of  the  tangent  component 


is  based  on  the  previous  projection  Vjjg  and  on  the  nor¬ 
mal  component  at  the  current  propagation  point.  Differ¬ 
entiating  this  equation  we  obtain 

3V, 


dV 


It 


'It  ^'^It  ^'^It 


=  -  cotci  dV.  +  -  da  (3.6) 

sina  OS  xn  ^.^2^ 

Note  that  the  error  in  depends  on  the  error  in  the  previous 
projection  (dV^g) ,  the  error  in  the  normal  component  at 
the  current  propagation  point  (dV^^) ,  and  the  error  in  the 
measurement  of  the  angle  a  (da). 

The  result  of  these  various  errors  is  that  when  the  pro¬ 
pagation  reaches  Aj^,  the  velocity  vector  attributed  to  Aj; 
by  the  propagation  procedure  will  differ  from  the  velocity  vec¬ 
tor  originally  computed  at  Aj^.  Therefore,  at  the  point  A^ 
we  compute  the  error  between  the  propagation  velocity  estimate 
Vj^  and  the  original  velocity  vector  Vj^  and  compute  the  error 


If  this  error  is  less  than  some  tolerance,  then  this  propagation 
procedure  is  stopped  at  point  otherwise  a  correction  proce¬ 
dure  is  applied.  If  we  consider  the  error  AV]^  as  having  been 
accumulated  in  the  previous  n  steps,  then  the  average  velocity 
error  in  one  step  is 


ibli 

so  we  have  ni*V^  as  the  velocity  error  at  the  m  step  and 
propagate  this  velocity  error  step  by  step  backward  to 
correct  the  estimated  velocity  vector  at  each  point  along 


the  same  contour. 


•  Experimental  results 

We  applied  the  corner  motion  estimation  and  velocity  pro¬ 
pagation  algorithms  to  two  sets  of  motion  pictures.  Section 
4.1  describes  the  corner  motion  estimation  results,  and  Sec¬ 
tion  4.2  describes  the  propagation  results. 

4.1  Corner  motion  estimation 

The  three  corner  motion  models  described  in  Section  2 
were  applied  to  two  image  sequences  containing  two  frames  each 
(Figures  6-7). 

We  first  describe  the  application  of  the  structural  model 
presented  in  Section  2.1.  Corners  are  "provisionally"  detected 
using  the  corner  detection  algorithm  described  in  Kitchen  and 
Rosenfeld  I  6 ] .  Next,  a  small  window  around  each  corner  is 
analyzed  to  obtain  a  more  accurate  description  of  the  corner. 
Based  on  the  assumption  that  the  corner  locally  contrasts  with 
its  surround,  a  local  thresholding  procedure  (Milgram  [  7  ] )  is 
used  to  segment  the  window.  The  corner  is  then  relocated  to  a 
maximum  curvature  boundary  point  in  the  threshclded  window. 

The  slopes  of  the  line  segments  meeting  at  the  corner  are  com¬ 
puted  using  a  one-dimensional  (slope)  Hough  transform  procedure 
(only  slope  need  be  computed  since  the  lines  are  constrained 
to  pass  through  the  corner  point. )  The  corners  detected  by 
this  procedure  are  marked  with  dark  crosses  in  Figures  6a  and  7a. 

To  overcome  the  effects  of  various  sources  of  error  on  the 
motion  estimation,  several  quadruples  of  line  segments  are  used 
to  compute  estimates  of  Ax'^,  Ay'^,  Ax'^  and  Ay'^,  with  the 
final  motion  estimate  taken  as  the  average. 


The  results  for  the  airplane  in  Figure  6  are  displayed  in 
Table  1.  The  estimated  motion  vectors  were  obtained  by  the 
authors'  examination  of  digital  enlargements  of  the  images. 

No  useful  results  were  obtained  for  the  moving  car  in  Figure  7. 
There  are  several  reasons  for  this: 

1.  The  grey  level  corners  in  the  car  are  much  more  rounded 
than  the  airplane^,  and  the  motion  estimates  are  sensitive 
to  the  corner  location;  and 

2.  The  spatial  resolution  of  Figure  7  is  not  high  enough  to 
allow  us  to  place  a  sufficiently  large  window  around  a 
corner  for  segmentation  which  does  not  contain  some  other 
image  feature. 

The  least-squares  corner  motion  estimates  presented  in  Sec¬ 
tion  2.2  require  that  we  first  compute  the  normal  component  of 
motion  along  the  contour  in  the  neighborhood  of  the  corner.  The 
magnitude  of  the  normal  component  and  the  components  of  the  unit 
normal  vector  on  the  x  and  y  axes  are 


follows : 


1  1  1 

I  =  —  Z  I  fl  -I  > 

X  4  k+«,,i+m,j+l  k+i,i+m,j 

.  1  ^  ^ 

Tssiryfl  -T  } 

y  4  j^_Q  k+«,,i+l,j+n  k+il,i,j+n^ 

I  =  -  Z  E  {I  .  -I  } 

t  4  jj^_Q  k+l,i+in,j+n  k,i+m,j+n 

The  unit  of  length  is  the  grid  spacing  interval  in  each  image 
frame  and  the  unit  of  time  is  the  image  frame  sampling  period. 

Tables  2  and  3  contain  the  results  of  applying  both  least- 
squares  corner  estimators  to  the  images  in  Figures  6  and  7, 
respectively. 
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8 
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1.0 
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Table  1.  Motion  vectors  for  corners  in  Figure  6. 
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1 
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-1.7 

1.0 
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1 

1 

-2.9 
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1 

1 

-2.4 

0.1 

-3.5 

3.7 

Table  2.  Motion  vectors  for  turning  points  of  Figure  6 
(airplane) 

(1)  Under  local  translation  assumption. 

(2)  Under  local  planar  motion  assumption. 
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41 
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21 

24 

-1.8 
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-1.4 

-0.5 
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39 
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-1.5 
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-1.4 

-2.5 

41 
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-1.3 
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-1.2 

-0.4 

3.0 

Table  3.  Motion  vectors  for  turning  points  of  Figure  7 
(traffic ) . 

(1)  Under  local  translation  assumption. 

(2)  Under  local  planar  motion  assumption. 


4.2  Velocity  propagation 

We  applied  the  propagation  technique  to  the  two  image 
sequences  displayed  in  Figures  6  and  7. 

The  propagation  technique  was  implemented  as  follows: 

1)  Velocity  vectors  are  first  determined  at  a  set  of 
"corner"  points  in  the  first  frame  by  the  least-square 
corner  motion  estimator  which  assumes  that  the  corner 
motion  is  a  2-D  translation. 

2)  The  velocity  vector  at  the  corner  is  propagated  along 
the  contours  that  meet  at  the  corner  until  a  second 
corner  point  is  encountered.  The  contours  are  followed 
by  a  very  simple  maximum  gradient  technique.  A  velo¬ 
city  vector  is  not  computed  at  every  pixel  on  the  con¬ 
tour,  but  only  at  every  pixel,  to  reduce  the  error 
in  a . 

3)  When  the  terminating  corner  point  is  reached,  the  propa¬ 
gation  is  stopped  and  the  error  velocity  vector  is  com¬ 
puted.  If  this  error  is  greater  than  a  preset  tolerance, 
then  the  error  velocity  vector  is  back-propagated  along 
the  same  contour. 


The  results  of  the  propagation  are  displayed  in  Figures  8-9. 
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Velocity  field  using  the  propagation  technique 
along  the  contours  of  the  moving  airplane  shovm 
in  Figure  6. 
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Velocity  field  using  the  propagation  technique 
along  the  contours  of  the  moving  car  shown  in 
Figure  7. 


5.  Summary 

We  have  presented  a  contour-based  approach  to  motion  estimation 
based  on  first  estimating  motion  at  image  corners  and  then  propa¬ 
gating  these  motion  estimates  along  image  contours.  One  potential 
advantage  of  such  an  approach  over  others  such  as  [3-4]  is  that 
motion  information  is  not  integrated  across  the  boundaries  of 
moving  objects,  but  only  along  such  boundaries.  Since  very  often 
the  only  reliable  source  of  motion  information  is  at  object  bounda¬ 
ries  (when,  for  example,  object  interiors  are  homogeneous)  it  is 
important  that  motion  estimation  techniques  yield  accurate  motion 
estimates  at  boundaries. 

The  examples  presented  in  Section  4  both  consisted  of  a  single 
object  moving  across  a  homogeneous  background.  The  propagation 
technique  presented  in  Section  3  would  need  to  be  modified  to  be 
applicable  to  more  complex  image  sequences  containing  multiple  moving 
objects  and  occlusions  so  that  motion  information  is  not  propagated 
from  one  object  to  another. 
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